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Abstract. In the radio regime the mass- loss rate of AGB stars is best probed 
using molecular (and atomic) line emission arising in the CSE formed by the stel- 
lar wind. The numerical modelling of the circumstellar emission where intricate 
interplays between physical and chemical processes take place, is a challenge. 
The derived mass-loss rates depend crucially on the assumptions in the cir- 
cumstellar model, of which some can be constrained if enough observational 
data exist. Therefore, a reliable mass-loss-rate determination for an individual 
star requires, in addition to a detailed radiative transfer analysis, good obser- 
vational constraints in the form of multi-line observations and radial brightness 
distributions. Of the methods used to estimate mass- loss rates from galactic 
AGB stars those based on radiative transfer modelling of CO line emission are 
most commonly used and possibly also the most accurate. Typically, CO multi- 
transitional observations can constrain the mass- loss rate to better than 50%, 
within the adopted circumstellar model. Comparison with complementary meth- 
ods, such as estimates based on dust radiative transfer modelling coupled with 
a dynamical model, are consistent within a factor of three. 



1. Introduction 

Mass loss is the single most important process during the final evolution of low- 
and intermediate- mass stars on the asymptotic giant branch (AGB). Its exis- 
tence and overall characteristics, e.g., magnitude, geometry and kinematics, are 
well established. However, much of its finer details, e.g., degree of asymmetry, 
non-homogeneity, and temporal variation are essentially unknown. This is un- 
fortunate since even a modest mass-loss-rate change by a factor of two will have 
a profound effect on the evolution of the star, its nucleosynthesis, and its re- 
turn to the interstellar medium. The mass loss creates a circumstellar envelope 
(CSE) of gas and dust around the star. The low temperature of the central star 
allows the formation of a wide variety of molecular species in its atmosphere, 
and the expanding gas is therefore mainly in molecular form. The chemistry in 
the CSE itself can be very rich and it depends on the C/O-ratio, the thickness 
of the envelope, and the strength of the ambient ultraviolet radiation field. 

Presented here is a short overview of various methods for estimating mass- 
loss rates based on observations in the radio regime. Special attention will be 
given to CO line observations which has been considered, when supplemented by 
detailed radiative transfer modelling, as a very (perhaps the most) reliable AGB 
mass-loss-rate estimator. Uncertainties involved in the modelling procedure will 
be discussed. 
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Figure 1. a) Mass-loss-rate distributions for the S-star (solid line; 40 stars; 
Ramstedt et al. 2006), M-star (dotted line; 77 stars; Olofsson et al. 2002; 
Gonzalez Delgado et al. 2003), and carbon star (dashed line; 61 stars; Schoier 
& Olofsson 2001) samples, b) Envelope gas expansion velocity distributions 
for the S-star (solid line), M-star (dotted line), and carbon star (dashed line) 
samples, c) Derived mass-loss rates plotted against the gas expansion ve- 
locities for the S-star (dots), M-star (squares), and carbon star (triangles) 
samples. 



2. Probing the mass-loss rate 



The discussion in this section will be limited to the more commonly used meth- 
ods in deriving mass loss rates for AGB stars based on observations of spectral 
line emission in the radio regime. Although, radio continuum observations could 
in principle be used to drive dust mass- loss rates, and through the use of a dy- 
namical model the gas mass-loss rates, they are mainly only included as part 
of more extensive multi-wavelength SED modelling which is usually better con- 
strained from observations in the infra-red (and optical) part of the spectrum. 
In Sect. 4.2 a comparison between gas mass- loss rates derived from CO line 
emission and dust continuum emission is provided. 



2.1. HI 21 cm line emission 

Hydrogen is the most abundant element in the material expelled from an AGB 
star. Molecular hydrogen (H2) is difficult to observe since it has no allowed 
electric dipole rotational transitions and its (weaker) quadrupolar transitions 
lie in a spectral range in the infra-red that is not accessible from the ground. 
The atomic hydrogen (HI) line at 21 cm could potentially be a useful probe 
of the gaseous envelopes around AGB stars. HI is also hard to ionize which 
would trace the CSE to large spatial scales. However, chemical models predict 
the gas in the wind to be atomic only when the central star has an eff ective 
temperature higher than about 2500 K (e.g. iGlassgold Hugginsl [l983). HI 



can also be produced from photodissociation of H2 by the interste llar radiation 
field at large radial distances from the star ( Gussie et al.1 [l995l ). There are 



some observational difficulties in that the emission is inherently weak and that 
HI emission from the interstellar medium provides an intense background from 
which the circumstellar contribution needs to be separated. 

Mass- loss rates based on HI 21 cm line observations are still very uncer- 
tain. This is partly because the HI abundance distribution may be compli- 
cated but there is also evidence that the AGB wind interacts with the ISM 
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at large radial scales. Current HI mass- loss rate estimates show a correlation 
with those obtained from other methods, such as CO line observations, but the 
discrepancy for an individual source can be as large as an order of magnitude 
dLe Bertre fc Gerard I200I : 1 Gerard Le Bertrdl2006h . 



2.2. CO line observations 

CO is usually the most abundant circumstellar molecule after H2 and multi- 
transition CO line observations, when supplemented by a detailed excitation 
analysis, are commonly regarded as one of the most accurate methods for deter- 
mining mass-loss rates on the AGB. The CO molecule is relatively easy to excite 
even in low temperature low density regions and, like H2, self-shileds against the 
interstellar uv-radiation field which combines to make it an excellent probe of 
the circumstellar medium. 

Several surveys of the mass-loss properties of AGB stars have been per- 
formed using CO line emission (e.g., Loup et al. 1993; Bieging & Latter 1994; 
Kahane & Jura 1994; Groenewegen et al. 1999; Hiriart & Kwan 2000; Schoier 
k Olofsson 2001; Olofsson et al. 2002; Ramstedt et al. 2006; Teyssier et al. 
2006). It is found that the derived mass-loss-rate distribution has a median 
value of ~ 2 x 10 M yr _1 with no obvious dependence on the photospheric 
C/O-ratio (see Fig. la). The mass-loss rates correlate well with the expansion 
velocity of the CSEs (Fig. lb) and there are some indication that carbon stars 
generally have larger expansion velocities than S-stars and M-type AGB stars 
(Fig. lc). The observed and derived properties of AGB CSEs and their corre- 
lation with stellar properties such as, e.g., luminosity and effective temperature 
lend support to the common concensus that these stellar winds are driven by 
radiation pressure on dust grains and that pulsation m ay play an important role 
(ISchoier & Olofssonl[200ll : iBereeat Sz Chevallied l2005h . 



Single-dish observations in the radio regime generally do not spatially re- 
solve the CSEs around AGB stars. The first systemat ic survey of AGB molecular 
envelopes at high spatial resolution was performed bv lNeri et al.l ( 19981 ) using CO 



J = 1 — > line emission using the Plateu de Bure interferometer. They concluded 
that the CSEs generally have an overall spherical symmetry and expand isotrop- 
ically at near constant expansion velocity. These results have further stimulated 
the continued use and development of the "standard model" described in Sect. 3. 

2.3. Maser emission 

Together with CO observations, maser line emission from molecules such as 
OH and H2O (and to a lesser extent SiO) has commonly been used to derive 
mass- loss rates for AGB stars. This procedure is biased toward AGB stars 
in which the photospheric C/O-ratio is < 1 which allows for the formation of 
high abundances of H2O and its photodissociation product OH. The excitation 
conditions requires the maser emission to emanate from a region relatively close 
to the central star allowing for detailed studies of the inner part of the CSE where 
the mass- loss is initiated and the wind accelerated. Maser lines are usually very 
brigh t and can be used to study AGB stars at larger (extragalactic) distances 
(e.g., Ivan Loon et al.|[200TI ) . 



Mass-loss rate estimates based on maser emission are inherently more un- 
certain than those obtained from CO measurements. Detailed radiative transfer 
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modelling indicates that the maser strengths are very sensitive to changes not 
only in the chemistry, density, temperature, and velocity structures of the gas 
but also in the exact composition of the dust particles pres ent in the wind (e.g., 
iGoldreich fc Scovilld [1979 : iBabkovskaia fc Poutaner3l2006h . 



3. The standard model 



The numerical modelling of CSEs, where intricate interplays between physical 
and chemical processes take place, is a challenge. This applies, in particular, 
to the important inner regions of the CSEs, where very currently few observa- 
tional constraints are available and our knowledge is very limited. Atomic and 
molecular line emission in the radio regime predominately probe the outer en- 
velope regions where many of these complexities can be neglected and a fairly 
simple, yet reasonably realistic, CSE can be assumed. Here a brief summary of 
the assumptions used in the "standard model" often adopted when estimating 
mass-loss rates from CO line observations in the radio regime. As an example a 
best fit model for the high mass-loss rate carbon star LP And is shown in Fig. 2. 
In Sect. 4 some of the uncertainties involved in this modelling procedure will be 
discussed. 



3.1. Basic assumptions 

The most commonly adopted starting point for detailed modelling is to assume 
the CSEs to be spherically symmetric, to be produced by a constant mass-loss 
rate, and to expand at a constant velocity. The density structure, pu 2 , (as a 
function of distance, r, from the central star) can then be derived from the 
conservation of mass 

M 

PR 2 = nu 2 m H2 = - — 5—, (1) 
4:7rr z v e 

where M is the hydrogen gas mass loss rate and v e is the gas expansion velocity 
taken from the observed line profiles. It is further assume d that the hydrogen is 
in mo lecular form in the region probed by the CO emission (Gl assgold fc Hugginsl 
1983). In the CO excitation analysis a CO/H2 abundance ratio needs to be 
adopted. Typically the assumed ratio is in the range 2 x 10~ 4 — 1 x 10~ 3 with 
M-type AGB st ars in the lower end and carbon stars in the upper end (e.g., 
ICherchnefj [20061 ). 

The CO molecules are generally not excited according to local thermody- 
namic equilibrium (LTE) in the regions of the envelope which contribute signif- 
icantly to the observed emission. The common approach is instead to assume 
that statistical equilibrium is prevailing. The excitation analysis then requires a 
detailed knowledge of collisional rate coefficients between CO and H2 molecules. 
Fortunately, CO-H2 is a relatively well studied system (see summary in Schoier et 
al. 2005b) and calculated collisional rate coefficients currently exist for temper- 
atures in the range 5 — 400 K and involve energy levels up to J = 29 and J = 20 
for collisions with para-H2 and ortho-H2, respectively. Extrapolations up to a 
collisional temperature of 2000 K involving energy le vels up to J = 40, adequate 
for AGB CSEs, are provided bv lSchoier et al.1 ( 2005bl ) and available for download 



though the on-line database LAMDA (www.strw.leidenuniv.nl/~moldata). 
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Figure 2. Multi-transition CO millimetre-wave line emission observations 
(histogram) of the carbon star LP And overlayed with model predicions (solid 
line). In the cooling panel the dotted and dashed lines are H 2 and CO line 
cooling, respectively. The solid line is adiabatic cooling due to expansion of 
the gas. In the temperature/optical depth panel the solid line is the excitation 
temperature of the CO J = 2 — > 1 transition and the dotted line is the kinetic 
gas temperature. The dashed line is the tangential optical depth of the CO 
J = 2^1 transition at the line center (figure from Schoier & Olofsson 2001). 



In AGB CSEs an ortho-to-para ratio of 3 is usually adopted when weighting to- 
gether collisional rate coefficients for CO in collisions with ortho- and para-HV 
Generally, and in particular for low mass-loss-rate objects, it is important that 
the excitation analysis also includes radiative excitation through the first vibra- 
tionally excited (v = 1) state at 4.6 ^m (see Sect. 4.1). The models also include 
both a central source of radiation and the cosmic microwave background radi- 
ation at 2.7 K. The central radiation emanates from the star itself, which may 
be approximated by a blackbody. In some sources, where a prominent infrared- 
excess is observed, fitting the spectral energy distribution with two blackbodies 
provide a first approximation to the radiation field that the CSE is subjected to. 
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Models that include an extend ed dust component in the CO excitation analysis 
may also be constructed (e.g.. ISchoier et al1l2002l ). 



The size of the CO envelope is an important parameter, and the derived 
mass loss rate will depend on this (Sect. 4.1). If no radial brightness distri- 
butions exist the size of the circumstellar CO envelope can be estimated using 
the modelling presented in Mamon et al. (1988). It includes photodissociation, 
taking into account the effects of dust-, self- and H2-shielding, and chemical 
exchange reactions. 

3.2. The energy balance equation 

Once the level populations are known the kinetic gas temperature can be cal- 
culated in a self-consistent way by solving the energy balance equation (e.g., 
Goldreich & Scoville 1976), 

where 7 is the adiabatic index, k the Boltzmann constant, H the total heating 
rate per unit volume, and C the total cooling rate per unit volume. The first 
term on the right hand side is the cooling due to the adiabatic expansion of the 
gas. Additional cooling is provided by mainly molecular line emission from CO, 
calculated from the derived level populations. In addition, HCN could be an 
important coolant in the inner parts of the envelope of carbon stars and H2O in 
the case of M-type AGB stars. 

The mechanism responsible for the observed mass loss is probably radiation 
pressure acting on small dust grains, which in turn are coupled to the gas. The 
radiation pressure on the dust grains will give them a drift velocity, v^ r , relative 
to the gas (Gilman 1972; Goldreich & Scoville 1976) As a result of the dust- 
gas drift, kinetic energy of the order of ^rri]i 2 v\ r will be transfered to the gas 
each time a particle collides with a dust grain. This is assumed to provide 
the dominating heating of the gas. Additional heating is provided by, e.g., the 
photoelectric effect (Huggins et al. 1988) that is important in the cool, tenuous, 
outer parts of CSEs around (in particular) high mass loss rate stars. 

When solving the energy balance equation several free parameters describ- 
ing the dust, i.e., via the dust-gas collisional heating, are introduce d. These are 

highly uncertain, but affect the derived line intensities. For example. ISchoier Qlofsson 



mg_ 
(20 



2001) define a new parameter that contains most of the these uncertainties, 

where $ is the dust-to-gas mass-loss-rate ratio, p<\ the density of a dust grain, 
ad the radius of a dust grain. The /i-parameter can be constrained if enough 
observational constraints are available (see Sect. 4.1). 



4. Error analysis 



A full, detailed, error analysis of the estimated mass loss rates is not possible due 
to the relatively large number of free parameters entering the model. Instead, 
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Figure 3. x 2 -maps illustrating the sensitivity of the standard model to the 
adjustable parameters; the mass-loss rate (M) and the /i-parameter. W Ori 
is a low mass-loss-rate object, R For is of intermediate mass-loss rate and LP 
And is a high- mass- loss rate object. The thick contour marks the 1-cr level of 
confidence (figure is adopted from Schoier & Olofsson 2001). 



some of the more important parameters may be varied in order to illustrate the 
sensitivity of the model, and to be able to get a rough estimate of the errors 
involved in the mass loss rate estimates. 



4.1. Sensitivity tests 

Detailed sensitivity tests within the adopt ed circumstella r model (see Sect. 3) 



have been performed by several authors (IKastnerl 119921 ; iGroenewegenl 11994 ; 
ISchoier fc Olofsson! [20011 : Olof sson et al. I l2002l ). In Fig. 3 chi-square contour 



plots for three example stars (W Ori, R For, and LP And), produced by varying 
the mass loss rate and the /i-parameter, illustrating the accuracy in the determi- 
nation of these two adjustable parameters when other parameters are held fixed 
(except for the CO outer envelope radius which is allowed to vary according to 
Mamon et al. 1988). 

LP And is a high mass-loss-rate object where the excitation of 12 CO is 
dominated by collisions. Consequently, the line intensities are very sensitive to 
the temperature structure, i.e., the /i-parameter. In the high mass loss rate 
regime the line intensities are also insensitive to the adopted mass loss rate 
since an increase in M leads to more cooling (CocM while ffocM ' 5 ), which 
compensates for the increase of molecular density. This "saturation" -effect of 
the line intensities for high mass loss rates has been noted before in other models 
where the cooling by CO is treated in a self-consistent manner (Sahai 1990; 
Kastner 1992). W Ori, is a low mass loss rate object with low optical depths 
in the lowest rotational transitions. In the low mass-loss-rate regime the line 
intensities scale roughly linearly with M. As expected the radiation emitted 
by the central star plays an important role in the excitation and the CO line 
intensities are less sensitive to the temperature structure (/i-parameter). R For 
presents an intermediate case between LP And and shares properties of both the 
low and the high mass-loss-rate stars. An intermediate mass- loss-rate model is 
generally more sensitive to the temperature structure than the radiation field, 
thus resembling the high mass loss rate objects. However, as in the low mass 
loss rate regime the line intensities scale roughly linearly with M. 



8 



Schdier 




10 6 10 5 10 4 



M [M yr ] (dust + dynamical model) 

Figure 4. Derived mass-loss rates from CO multi-line modelling and from 
SED modelling using a dust radiative transfer code coupled with a dynamical 
model. The sample consists of 5 M-type AGB stars (squares) and 5 carbon 
stars (triangles). The solid line marks the 1:1 correlation and the dotted lines 
a factor of three deviation (Ramstedt et al. 2007, in prep.). 



In additional tests, ISchoier &: Olofssonl ( 200ll ) noted that the assumed en- 
velope size r p does not significanly affect the derived mass loss rate in the high 
mass loss rate regime since the density is too low to excite the CO molecules 
effectively (the lines are sub-thermally excited) in the cool outer parts of the 
envelope (cf., Fig. 2), i.e., the emission is excitation limited. In contrast, for 
low mass-loss-rate objects the size of the envelope is important when determin- 
ing the mass loss rate, i.e., the emission is pho todissociation limited (at least 
for the lowest transitions). ISchoier &: Olofssonl (|200ll ) made tests which show 
that for typical envelope parameters the transition from photodissociation- to 
excitation-limited emission lies at ~5 xlO _7 M0yr _1 for the J=l—>0 transition. 
Higher J-lines sample hotter and denser gas, located closer to the star, and are 
therefore less sensitive to the choice of the envelope size. 

In summary, it is concluded that within the adopted circumstellar model 
the estimated mass loss rates are accurate to about ±50% (neglecting errors 
introduced by the uncertain CO abundance and the distance estimates) when 
good observational contraints are available, but one should keep in mind that 
the causes of the uncertainty varies with the mass loss rate. 

4.2. Comparison between various authors and methods 

A rough way to estimate the errors involv ed is to compare w i th th e results 
obtain ed from other self-co nsistent models. Schoier Sz Olofsson ( 200ll ) and re- 



cently [^£SsieL^t^ah ( 20061 ) present comparisons of derived mass-loss rates from 



detailed modelling found in the literature for a large number of sources. They 
usually find a very good agreement, within 30%, when adjustments for differ- 
ences in adopted CO abundance and distance have been made. For an individual 
star, however, the discrepancy can be as large as a factor of three. 

Dust radiative transfer modelling when coupled with a dynamical model 
of the wind can provide an independent estimate, to that of CO, of the gas 
mass-loss rate. In Fig. 4 derived gas mass-loss rates obtained from these two 
complementary methods are compared for a sample of intermediate- to high 
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Figure 5. left — Interferometric observations of CO line emission towards 
the M-type AGB star R Cas obtained by the PdBI (Castro-Carrizo et al. 
2007, in prep.), right - The azimuthally averaged brightness distribution of 
the CO line emission around R Cas and the derived temporal modulation of 
the mass loss rate (Schoier et al. 2007, in prep.). 



mass-loss-rate objects (Ramstedt et al. 2007, in prep.). It is found that the two 
methods are generally consistent within a factor of three. Possible systematic 
effects here are the adopted CO fractional abundance and dust properties. 

4.3. Departures from the standard model 

The sucessful modelling of the vast majority of AGB stars, assuming a spher- 
ically symmetric CSE, suggests that pronounced axi- or non-symmetric mass 
loss is not a common phenomenon among AGB stars. However, there are cases 
wher e clear bipolar outfl ows have been detected such as the carbon star V Hya 
(e.g., iHirano et all 120041 ). 



A small fraction (~ 10%) of all carbon stars have detached molecular shells 
produced during a relatively short period of inte nse mass loss thought to be re- 
lated to thermal pulses (e.g.. lQlofsson et "ai]|l996l ). To be able to put constraints 



on the mass loss and its variation with time, high spatial resolution observa- 
tions are needed. Radial brightness distributions also reveal smaller modula- 
tions of the mass- loss rate for a number of AGB stars (e.g., Schoier & Olofsson 
2001; Teyssier et al. 2006; see also Fig. 5). Moreover, an interacting wind sce- 
nario where the AGB wind sweeps up the ISM need s to be taken into account 



when interpreting m ass-loss rate modulations (e.g.. iLe Bertre Sz Gerard! 12004k 



Schoier et al J 12005a! ) and the degree of dumpiness of the circumstellar medium 



needs to be fully investigated. 



5. Conclusions 

The mass loss rate determination for AGB stars depends crucially on a number 
of assumptions in the CSE model. A reliable mass loss rate determination re- 
quires, in addition to a detailed radiative transfer analysis, good observational 
constraints in the form of multi-transition observations and radial brightness dis- 
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tributions. A standard CSE model, assuming a single smooth expanding wind 
produced by a continuous mass loss, can explain multi-transition single-dish line 
observations in the majority of AGB stars and objects with clear deviations 
from the basic model assumptions can easily be identified. However, in order 
to investigate small deviations in terms of mass-loss rate modulations (typically 
lower than a factor of two), departure from spherical symmetry, or the degree 
of dumpiness, high-spatial-resolution interferometric observations are required. 
Such new observations are already becoming available and reveal the compo- 
sition of AGB circumstellar envelopes are usually pathcy and clumpy to some 
degree also at larger radial scales from the central star. The next generation of 
radiative transfer codes, as well as dynamical and chemical models of circum- 
stellar envelopes, need to take such effects into account. 

Acknowledgments. F.L.S. acknowledges financial support from the Swedish 
research council. 
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